Mon. Not. R. Astron. Soc. 000, 0-|l2| (1996) Printed 1 February 2008 (MN M£X style file vl.3) 



A homologous recycling model for hot galactic coronae 



Alexei G. Kritsuk 1 ' 2 ' 3 * 

1 Institute of Astronomy, University of St. Petersburg, Stary Peterhof St. Petersburg 198904, Russia 

2 Max- Planck- Institut fur Astrophysik, Postfach 1523, D-85740 Garching, Germany 

3 Max- Planck- Institut fur extraterrestrische Physik, Postfach 1603, D-85740 Garching, Germany 



Accepted 1995 December 11. Received 1995 July 6; in original form 1995 January 3 



On 
On 

(X, 
(N 



> 

in 
o 

(N 
O 

VO 
On 

Oh 
6 



13 



ABSTRACT 

An equilibrium model is presented for a hydrostatic isothermal hot gas distribution in 
a gravitational well of a giant elliptical galaxy immersed in a massive dark halo. The 
self-consistently determined gravitational potential of the system provides the optical 
surface brightness distribution of the galaxy, matching the de Vaucouleurs i? 1 / 4 -law, 
if the sound speed in the gas c, stellar velocity dispersion <7„, and velocity dispersion 
of dark matter ctdm are related by <tdm — c ~ v2<7*. While the first equality follows 
from virial considerations, the second one relies on the observed similarity of optical 
and X-ray brightness profiles. The thermal equilibrium of the gas is described in the 
framework of a one-zone model, which incorporates radiative cooling, stellar mass loss, 
supernova heating, and mass sink due to local condensational instability. The sequence 
of saddle-node and saddle- connection bifurcations provides a first-order phase transi- 
tion, which separates the stable hot phase in the cooling medium for a reasonably high 
efficiency of condensation and precludes catastrophic cooling of the gas. The bifurca- 
tions occur naturally, due to the shape of the radiative cooling function, and this holds 
for a wide range of gas metallicities. The one-zone model implies a scaling relation for 
the equilibrium stellar and gas densities p* oc p 2 , which allows a hydrostatic, thermally 
stable distribution for the hot isothermal coronal gas. The resulting homologous re- 
cycling model reproduces the basic optical and X-ray properties of relatively isolated 
giant elliptical galaxies with hot haloes. Implications for cD galaxies in X-ray dominant 
clusters and groups are briefly discussed in the general context of galaxy formation and 
evolution. 
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1 INTRODUCTION 

The analysis of the extended X-ray emission of early-type 
galaxies has led to an extension of the cooling flow mode l 
[see Fabian, Nulsen & Canizares (1984) and Fabian ( 1994 )1 
applied to the class of bright elliptical galaxies (Thomas 
et al. 1986). The typical X-ray luminosities of such galaxies 
are in the range 10 39 — 10 42 erg s _1 , the gas temperature 
is 5 x 10 6 — 2 x 10 7 K, and the electron number density 
in the centre is 0.01 — 0.1 cm -3 . The galaxies usually con- 
t ain 10 9 - 1O 1O M of the gas in their extended hot haloes 
(Forman, Jones & Tucker 1985). Observations reveal a tight 



ilarity of X-ray and optical light distributions (Trinchieri, 
Fabbiano & Canizares 1986; Killeen & Bicknell 1988). 

The source of the gas in hot haloes of early-type galax- 
ies is usually attributed to the normal stellar mass loss with 



a specific rate p*/p* 



-2xl0~ 12 yr~ 



Gas injection from 



connection of X-ray and optical characteristics of the galax- 
ies, including a strong correlation between the X-ray and 
optical luminosities, al t hough with a conside rable scatter 
(Canizares et al. 1987; Donnelly et al. 1990), and a sim- 
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stars and su pernova events in the galaxy make the gas inho- 
mogeneous ( Mathews 199C ) and thermal instability causes 
matter to co ol out at radii r < r coo i, where t coo i < Hq 1 
( |Nulscn 198e| ). The dynamics of the gas, therefore, is deter- 
mined by the relative rates of the processes, controlling its 
mass and energy balance. Theoretical work done so far has 
demonstrated a variety of possible behaviours for the gas, 
ranging from supe rsonic galactic winds to su bsonic inflows 
[see Thomas ( jl986|), Ha ttori, Habe & Ikeuchi ( |l987| ), David, 
Forman & Jones (199C), Ciotti et al. (1991) and references 
therein]. However, due to the non-linear nature of radia- 
tive cooling, numerical modelling of inflows usually ends up 
with a cooling catastrophe in the centre of the galaxy, which 
overproduces the observed X-rays and precludes the study 
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of further evolution. As there are no certain observational 
indications of such a catastrophe, it is important to under- 
stand the nature of possible physical damping mechanisms, 
which are able to prevent it. 

It is an objective of this paper to present a purely hy- 
drostatic, thermodynamically stable model for gas recycling 
in an elliptical galaxy, which satisfies the observational con- 
straints. The model incorporates all the same physical pro- 
cesses as the 'standard' multiphase cooling flow model, and 
may also be useful in interpreting the X-ray emission from 
the central 'cooling flow' regions in dynamically quiet X-ray 
dominant clusters or groups, where the outskirts of the cen- 
tral dominant optical galaxy extend to the cooling radius. 



2 UNDERLYING GALAXY MODEL 

Let M(r) be a spherically symmetric configuration of the 
gravitating mass, including the stellar component of a 
galaxy, i Lire hut gas, and Lhe dark halu. Then Lire hydiu- 
static g aa diotribution in the gravitational potential 4>( r ) of 

M(r) can be described by 



kT /d hip d InT 



/d In 
V~dr 



GM(r) 



(1) 



p,m p \ d r d r / 

For the isothermal gas this reduces to 
kT d In p _ , 

j- - = -V0. 2 

pm p a r 

Here T and p are the gas temperature and density, p the 
mean molecular weight, m p the proton mass, and k Boltz- 
mann's constant. 

The stellar density distribution of the galaxy, p*, re- 
sponding to the same gravitational potential, satisfies 



2 [ d lnp» d lna r 2r] 

a T — - 1 1 = -V0, 

d r d r r 



(3) 



where <r r is the radial velocity dispersion of stars in the 
galaxy, the anisotropy parameter r\ is defined as r\ — 1 — 
of /of , and the subscript 't' represents the tangential veloc- 
ity component. In the case of an isotropic and uniform stellar 
velocity distribution (at — a r = const) it simply follows that 



id In p* 
d r 



-V0, 



(4) 



where cr, is the one-dimensional velocity dispersion. For an 
isothermal gas and an isotropic, constant stellar velocity dis- 
persion the stellar and gas densities are related as 



d In p 
d Inp* 



Cp'l 



pm p a t 
kT 



(•») 



(6) 



where c is the isothermal sound speed and C is the inte- 
gration constant. This relatio n provides the basis for the 
so-ca lled isothermal /3-model (Cavaliere & Fusco-Femianc 



197f), extensively used to describe the X-ray surface bright- 



ness of the hot gas in clusters of galaxies. 

Given equation (^) and assuming similar profiles for the 
gas density p and for the dark matter density pdm (and, 



therefore, ctdm = c)[j], one can define the total mass distri- 
bution self-consistently, using Poisson's equation 



--(r 

r 2 dr V 



i d In pt 
dr 



4ttC?. . 
— (p» + p + Pdm ) • 



(7) 



When j3 — 1/2 and in the centre p* : o is somewhat higher 
than pdm,o, the stellar system has a density distribution 
which is similar to the profile of a self-gravitating isothermal 
sphere, but has a density cut-off at the core radius of the 
dark matter halo, where the dark m atter density begins to 



fall as pdm oc r , cf. Burkert (1994). It turns out that such a 



solution of (|7j) gives the projected stellar mass density, which 
follows the de Vaucouleurs 7? 1//4 -law over several orders of 
magnitude in surface brightness, in a wide range of radii 
_0AR c < r < 1.5-Ro, where R c is the effective radius (Burkert 



1994). This can be expected if the stellar system formed in 
a dissipative process from post-infall shock heated gas with 
the star formation rate oc p 2 . 

The volume density of the j3 = 1/2 stellar syst em falls 



more steeply than p* oc r outside the core region ( Burkert 



1994), therefore the restriction on stellar velocity dispersion, 
imposed to derive ™, can be formulated as 



d lncr* 



d lnr 



+ V 



< 1, 



while the gas must be isothermal up to 



d InT 



d lnr 



<C 1. 



(8) 



(9) 



Due to projection effects it is difficult to demonstrate 
whether the above conditions are indeed satisfied in observed 
systems or not. However, in most cases early-type galaxies 
do exhibit rather flat line-of-sight velocity dispersion profiles 
and observations also show that the gas is quite isothermal. 
The details of observational data will be briefly discussed in 
Section |^. 



3 ONE-ZONE MODEL REVISITED 

A simple semi-analytical one-zone model allows an insight 
into the non-linear thermodynamics of the hot gas in the 
halo of a galaxy. While t he stability of gas equilibria was 
studied in Kritsuk (1992) with an emphasis on limit cycle 
solutions, considerations here concentrate mainly on classi- 
fication of bifurcations, which lead to the separation of a 
stable hot phase of the thermally unstable cooling medium. 
In the following subsections t he ba sic formulae of the model 
are compiled from Kritsuk (1992) in order to lead up to 



the bifurcation analysis. The only new physical process in- 
troduced is the feedback heating. It was neglected before, 
but may play an important role in highly non-equilibrium 
situations. 



t The approximate equality of the dark matter velocity disper- 
sion and the sound velocity of the gas has been demonstrated, 
e.g. by the simulations of X-ray clusters of galaxies in the cold 
dark matter cosmogony for those clusters that are not merging 
objects, experi encing a transient boost in the velocity dispersion 
of the system ( Navarro, Frenk fc White 1995| ). 
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3.1 Input physics 

The set of physical processes incorporated in cooling flow 
models includes radiative cooling, mass deposition due to 
thermal instabilities, supernova (SN) type la heating, and 
stellar mass loss for the galaxy, which is immersed in a dark 
halo. Here only the basic formulae are presented; one can 
find mo re ex tensive discussion elsewhere [e.g. in Sarazin & 
White (1987), see also references below]. The standard no- 
tation is kept where possible. 



3.1.1 Stellar mass loss 

The efficiency of gas deposition due to stellar mass loss and 
SN la explosions is proportional to the stellar mass density: 
ap*(r), where a — a, + a sn = const incorporates contribu- 
tions fromj2Mre_cjuj^ 

s _1 ( |Faber fc Gallagher 197(1 |Sarazin fc White 1987| ), and 
mass loss by supernovae ct sn — 1.33 x 10~ s~ (the esti- 



mate corresponds to r sn = 0.24 h 50 SNu, M sn = 1.4 Mq, 
and (M/Lb)* = 8M0/L0). The supernova type la rate, 
r sn , is uncertain by a factor of the or der of 1.5 according to 
van den Bergh & Tammann (1991), where /2.50 = Ho/50 
km s~ 1 Mpc _1 , but a significantly lower estimate r sn = 
0.06+ 0.03/^0 SNu is given by Turatto, Cappellaro & Benetti 
(|1994[). 



3.1.2 Thermal instabilities 

Drops of condensed material cool out from the hot gas be- 
cause of thermal instabilities at the rate of p« = b\(n)p, 
where the dimensionless parameter b = const £ [0, 1] con- 
trols the efficiency of condensation. The Heaviside function 
X and the instability growth rate n are defined as 

n, if n > 0: 



X(n) 



0, otherwise; 



d ( ln 5 p\ _ 1 f 2 P A ^dA^ ap. 



dt 



\ T 



AT) 



(10) 



(11) 



Here c v is the specific heat at constant pressure, the first 
term in the rhs of equation ( |H| ) coincides with Field's insta- 
bility criterion, and the second d escribes stabilization due 
to stellar mass loss ( [Kritsuk lgggh . Note that if p* oc p 2 as 
in the (3 = 1/2 case, then p« oc p , which naturally implies 
a star formation rate from the cooled material proportional 
to the gas density squared. 



3.1.3 Radiative cooling 

Non-equilibrium radiative cooling as a function of tempera- 
ture and metallicity, A(T, Z) erg cm 3 g -2 s -1 , is c ompiled 



from the calculations of Sutherland & Dopita (1993), where 
s olar abundance ratios were taken from Anders & Grevesse 
(1989). Effects of dust on the cooling are not taken into 



account here since the coronal phase is believed to be es- 
sentially dust-free. If the gas-to-dust ratio were similar to 
Galactic values, and the dust were well mixed in the hot in- 
terstellar medium, it could be the dominant cooling ag ent at 
gas temperatures above 2 x 10 6 K (|Burke fe Silk 1974). This 



would essentially change the shape of the cooling function 
and therefore the gas stability properties outlined below. A 



crude estimate shows that dust cooling can be safely ignored 
if the dust content is ^ 0.01 of the Galactic value. 



3.1.4 Heating 

The rate of heating due to thermalization of stellar winds 
and due to type la SN events is assumed to be equal to 
ap,To, where To is the characteristic temperature of the heat 
source, To = (a*T„ + a sn T sn )/a. Heating by stellar motions 
is described by T* = um p o 2 /k = 6.47 x 10 6 K (an estimate 
is made for o, = 300 km s _1 ), and SN temperature T En = 



Jim,, v ■ 



1.09 x 10 9 K for E sn = 6 x 10 50 erg, M S1 



1.4 



3fc 

M Q . With these values T ~ 3.6 x 10' K, while for r sn = 0.06 
SNu the source temperature is lower: To ~ 1.7 x 10 7 K. 



3.1.5 Feedback heating 

Energy feedback by type II (lb) SN explosions due to 
star formation in cooled material can serve as an addi- 
tional source of gas heating. This slightly modifies formu- 
lae for a and To: a = a, + a sn + aib, where aa> = Etbpti, 
To = (a»T. + QfenTen + a fb T fb )/a. Here T fb = = 

3.6 x 10 9 K, Eth = 0.72E in the Sedov stage, and the su- 
pernova explosion energy E = 2 x 10 jl e rg. T he optimistic 
estimate £fb = 0.01 is given in Silk et al. (1986), assuming a 
rate of type lb SNe of 1 supernova yr -1 per 100 Mq yr~\ 
and the mass of ejecta returned to the ISM M c j = 1 M©. 
The associated mass feedback is negligible, i.e. a 3> EfbQ. 



3.2 Basic equations 

The mass and energy balance of the hot gas in the halo of 
a galaxy may be described locally by the autonomous set of 
ordinary differential equations: 

dp 



= ap» - pti 



(12) 
(13) 



dt 

— = - [«p,(T - T) - ( 7 - l)p ti T - p 2 A/c v ] 

where t is the time variable and a, p*, and To are treated 
as constants; c v is the specific heat at constant volume. The 
evolutionary solutions of equations ( |l2| ) and ([[3]) may be 
represented as trajectories in the (p,T) phase plane. The 
steady states of the system correspond to the fixed points 
and may be readily found as 



P = PeqCQ 



T = OoqTo, 



ap, 



(7 - l)Ao ' 



(14) 
(15) 

where Co = c(To) is the isothermal sound speed, Ao = 
A(To,Z), metallicity Z is assumed to be constant, and p eq 
and O cq are dimensionless functions of b, To, and Z: 



78 



1 l)(2-i/ eq ) = l + i 



(16) 
(17) 



Pcq = (1 - 7©eq)/Aeq. 

The subscript 'eq' denotes an equilibrium state, v = 
dlnA/dlnT, and A = A/A . 

The characteristic time-scale for evolution near the 
equilibrium, 
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(18) 



is a geometric mean of the local cooling time t CQ = 
c o/[(7 ~ l)pAo] an( l the time-scale for stellar mass loss 
t a = {ap*/p)~ l . 

For a given metallicity of the gas Z, stellar mass loss 
rate a, temperature of the heat source To, and efficiency 
of condensation b, one can find the equilibrium solutions, 
if the cooling function is known. It is a straightforward 
piece of algebra then to classify these equilibria accord- 
ing to their stability properties, using the standard linear 
technique (see Kritsuk 1992). This analysis relies on poorly 
known high-order derivatives of A with respect to tempera- 
ture and, therefore, the fine details of stability changes can- 
not be well resolved. At the same time, the main stability 
features, controlled by the nature of radiative cooling (res- 
onance line emission, producing a prominent maximum in 
the range of temperatures 10 4 ' 2 K ^ T Ss 10 7 ' 2 K, and ther- 
mal bremsstrahlung for T J> 10 7 ' 2 K), do provide physically 
important effects, which are discussed in the next section. 



3.3 First-order phase transition in a system with 
no feedback 

The efficiency of condensation b is the only essentially free 
parameter of the model. It is important to understand how 
the equilibria and their stability follow the changes of this 
control parameter. 

This is illustrated with a plot of the equilibrium tem- 
perature T eq against b for a fixed value of the heat source 
temperature To = 3x 10 7 K and solar metallicity (see Fig. |l|). 
There are no steady-state solutions when the condensation 
efficiency is low, b < b\ — 0.0018. The catastrophic cool- 
ing is unavoidable even in a system with reasonably strong 
SN heating. At 6 = bi, as a result of a fold bifurcation, a 
saddle-node pair of fixed points appears simultaneously, pro- 
viding stable and unstable steady-state solutions, classified 
as a node (higher temperature state) and a saddle (lower 
temperature state), respectively. The fold bifurcation is the 
simplest form of local bifurcation in that it requires only 
a one-dimensional phase space, only one control parame- 
ter, and only lowest-order terms (see, e.g., Guckenheimer & 
Holmes 1983; Thompson & Stewart 1989). It is structurally 
stable and, therefore, the small perturbation of the dynam- 
ics does not alter the bifurcation qualitatively. Equation ( |l6| ) 
shows that the condensation efficiency b is a unique function 
of the equilibrium temperature oq and depends solely on 
the logarithmic temperature derivative of the cooling func- 
tion f e q. Therefore the fold cafastrop/i^] occurs in the model 
due to the local structure of A(T). As the control parameter 
becomes larger, multiple steady states appear and annihi- 
late due to a series of bifurcations of the same kind. These 
further local changes in the phase portrait of the system at 



b\ < b < 0.054 are not so important because the new sta- 
ble solutions have negligibly small domains of attraction (or 
basins) in the phase plane. The most essential event in the 
range of b £ [0, 0.6] happens at b — bi, where the first pair of 
fixed points spontaneously appears. Although these changes 
in the dynamics do not preclude the catastrophic cooling, 
the birth of the saddle leads up to further transformations 
in the system. 

The next important event takes place at b 2 = 0.092 
due to a spontaneous birth of a stable limit cycle. This is a 
homoclinic connection bifur cation or a blue sky catastrophe 
for a periodic limit cycle (|Guckenheimer fc Holmes 198^ 



Thompson & Stewart 1989). The blue sky catastrophe is a 
global bifurcation^] in the phase plane in which the attract- 
ing limit cycle disappears (as b decreases) by collision with 
a saddle. The bifurcation is also structurally stable, that is, 
the full control-phase space diagram may be deformed, but 
qualitatively unchanged by small perturbations of the gov- 
erning equations. At b > b 2 the basin of the attracting limit 
cycle grows, while the amplitude of non-linear temperature 
oscillations becomes smaller. As a result, at b > 0.2 only 
small-amplitude stable oscillations survive, which disappear 
at b — 63 = 0.427 due to a supercritical Hopf bifurcation. 
Fig. |2] shows the huge basin of the stable spiral point at 
b — 0.5. This means that, whatever the initial temperature 
and density of the gas are, the thermodynamic state of the 
gas settles to an equilibrium with temperature ~ 10 7 K (for 
To ~ 10 7 - 5 K) in a finite time of the order of t s (see equation 
|is|). What actually happens at the threshold value 62 is a 
first-order phase transition, which separates the hot phase 
in the cooling thermally unstable medium. Thus, one has a 
stable thermodynamic description for a gas with embedded 
cold condensations for b > & cr it ~ 0.2, when the oscillations 
are small. The occurrence of the relaxation oscillations in the 
vicinity of a saddle-connect ion fo llows from theorem 6.1.1 in 
Guckenheimer & Holmes (1983). A detailed discussion of 
the cooling function properties, whic h res ult in the saddle 
connection, can be found in Kritsuk (1992). 

The key equation ([Hi]) does not contain an absolute 
value of the cooling, therefore one can conclude that, al- 
though the cooling is less efficient for metallicities lower than 
solar, the structure of equilibrium states, and, in particular, 
the equilibrium temperature T(b, Z) does not change dra- 
matically with Z. This is actually the case, as follows from 
the calculations based on the zero-field non-eq uili brium ion- 
ization cooling functions of Sutherland & Dopita (1993). For 
the same value of To as before and Z = 0.1 Zq the critical 
points, corresponding to the same catastrophic transitions, 
are 61 = 0.001 and b 2 = 0.097. The Hopf bifurcation oc- 
curs somewhat earlier, at 63 = 0.136. In the zero metallicity 
case 61 = 0.000 48 and b 2 « 0.029, but the nature of the 
saddle connection at 62 is different [positive trace case, see 



t The notion of catastrophe is introduced in non-linear dynamics 
to specify a spontaneous birth of an equilibrium state at a critical 
point. A fold is an example of such a discontinuous bifurcation. Its 
normal form, which embodies the local low-order approximation, 



may be written as x 



where for control |i > a pair 



of equilibria exists, one attracting and the other repelling. For 
u < the sol utions diverge to x — > — 00 as t — * 00 (Thompson & 
Stewart 19891). 



8 The profound changes in qualitative behaviour of the system 
near the threshold value b 2 are associated solely with the global 
topological configuration of the insets and outsets (the incom- 
ing and outcoming critical trajectories) of the saddle fixed point, 
which requires a finite region in the phase plane and cannot be 
described locally. A saddle connection bifurcation is known to oc- 
cur, for example, in the averaged autonomous van der Po l system 
[see equation (2.1.14) of Guckenheimer & Holmes (1983)1. 
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0.0001 0.0010 0.0100 0.1000 1.0000 10.0000 

Efficiency of Condensation, b 

Figure 1. Bifurcation diagram for To = 3 X 10 7 K and Z = Zq. The solid line shows the equilibrium temperature T = C qTo as a 
function of condensation efficiency b. Dotted lines indicate minimum and maximum temperature values for limit cycle solutions. The 
phase transition occurs at b = b^ • 




0.01 0.1 1.0 10.0 

Density 



Figure 2. Phase portrait for one-zone model with Tq = 3 X 10 K, b = 0.5, and Z = Zq. Solid lines with arrows show trajectories in the 
phase plane (temperature and density are dimcnsionlcss). Dashed lines indicate the insets of the saddle or separatrices, which separate 
the attraction domain of the stable spiral point. 
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Guckenheimer & Holmes ( 198E ) for details], and therefore 
an unstable limit cycle appears before the second key bifur- 
cation, at b < 62. These details, however, do not change the 
general picture of the phase transition due to the structural 
stability of bifurcations involved. 

The main conclusions of this subsection are the follow- 
ing: (i) the distributed sources and sinks of mass and en- 
ergy in the hot gaseous halo may be locally balanced and 
the thermodynamic equilibrium is stable if b > 6 C rit(To, Z); 

(ii) non- equilibrium solutions settle to the steady state from 
any initial conditions in its huge domain of attraction, and 

(iii) the equilibrium gas density is proportional to yfpl. It 
is worth emphasizing here that the phenomenon of phase 
transition in the system is not due to SN heating, but is 
due to a certain combination of the heating and distributed 
mass deposition. 

3.4 The role of feedback 

As shown in the previous subsections, the stable equilibrium 
can be reached by the system with no feedback heating. If 
star formation occurs in the condensed material, a relatively 
small energy feedback can modify the dynamics. This prob- 
lem will be studied in detail elsewhere. Here we concentrate 
only on the effects that feedback heating has on the positions 
of the fixed points in the phase plane. It is assumed that the 
feedback does not alter the growth rate of condensational 
instability (jll]). Then equation ( |l3| ) must be rewritten as 

dT 1 

— = -{ap f (T -T)+ p ti [Ti - (7 - m - P 2 A/cv} , (19) 

to account for the additional heating. The new quantity 
Ti — £ft,Tfb is introduced as the feedback temperature. It 
is easy to check that near the equilibrium equation (|l9]) is 
equivalent to (|l^) with To +Ti substituted for To in equation 
(tL3T) - Therefore the presence of additional feedback heat- 
ing does not change the structure of equilibria; only the 
characteristic heat source temperature increases, providing 
a higher equilibrium temperature for the gas. The role of 
feedback is somewhat more complex in situations that are 
far from equilibrium due to a time delay of the feedback star 
formation. 



4 HYDROSTATIC RECYCLING OF 
ISOTHERMAL GAS 

Perhaps 90 per cent of the papers on cooling flows contain 
an explanation for the flow's existence as due to a short cool- 
ing time-scale in the central region of a cluster or a group, 
where a cD or giant elliptical galaxy is located. It will be 
shown in this section that, while this condition is neces- 
sary for the flow to occur, it is not always sufficient. Note 
that the model discussed here is based on the same physics 
that is incorporated in most of the cooling flow models. The 
only difference is that special attention is paid to make the 
description more self-consistent and to incorporate the sta- 
bility analysis in order to restrict the parameter space and 
to avoid problems in numerical simulations^]. 



4.1 Does cooling gas necessarily flow? 

Simple comparison of equations (^|) and (^|) demonstrates 
that a self-consistent hydrostatic model for isothermal gas 
recycling in the hot corona of a giant elliptical galaxy can 
be constructed if 



P(<t» ,T ,Z, b) 



(20) 



i.e. in the same case in which the p„ distribution provides 
the observed optical surface brightness profile (see Section 
^|). The gas metallicity, Z, must be uniform over the galaxy. 
Equation (^o|) can be used as an additional constraint to 
reduce the dimension of the parameter space of the model. 
For example, the definition for /^(equation |B|), equation (|l5[), 
and the equilibrium condition (hq) allow one to calculate the 
value of b for given To , a,, and Z as 



kT 



27/im p (T» 



1 (2 



Finally, the integration constant C in equation 
obtained from equation (|l4|): 



C = 



■ 7c 4 



( 7 -l)A(T,Z)' 



(21) 



can be 



(22) 



The rate of recycling is a strong function of radius and 
scales as t^ 1 ~ pV 2 ~ p in the equilibrium. As the star for- 
mation rate in the cooling gas is proportional to the stellar 
density in the galaxy, the newly born stars and the associ- 
ated heating are distributed in such a way that they do not 
distort the structure of existing gaseous and stellar compo- 
nents. Therefore the recycling is homologous. 

For a comparison with cooling flow models one can for- 
mally calculate the mass deposition (accretion) rate for the 
recycling model as 



M(r) = / 47rr 2 p t i(r)dr = aM t {r). 



(23) 



For r 3> r c one has M ~ r , if p* ~ r~ n , and, thus the 
model predicts a scaling M oc r, if the stellar density dis- 
tribution is that of the isothermal sphere. Generic slopes of 
optical surface brightness profiles for giant elliptical galaxie s 
are in the range from 1.65 to 2.1 (Mihalas & Binney 1981). 
Therefore -0.1 & dlogM(r)/dlogr <, 0.4. While giant el- 



lipticals have relatively steep surface brightness profiles and 
correspondingly flat M(r), cDs, having flatter brightness dis- 
tributions, must have steeper mass deposition profiles. 

If the radiative cooling function A(T, Z) is assumed to 
be known, the list of input parameters of the model includes 
the following values: Pdm , P* > a *> a > To, and Z. The model 
could be applied in several consecutive steps. First, the ratio 
(Pdm + P )/P*o an d P*a must be determined as the best- 
fitting values for the Poisson's equation solutions to adjust 
the optical surface brightness profile. Secondly, the temper- 
ature and metallicity, derived from the spectral analysis of 
the X-ray data, provide the condensation efficiency b for the 
assumed value of To. An iterative procedure must be applied 
to get the value of To that satisfies the observed X-ray lumi- 
nosity. If the iterations converge, finally, all involved density 



^ For example, if non-linear relaxation oscillations of the ther- 
modynamic state of the gas occur, special care must be taken to 



ensure the convergence of the numerical method in use, cf. Friaca 
(|l993|). 
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distributions are known and the projected distribution of the 
X-ray emission may be compared with the observed one. 

Evidently, it is not always possible to find a hydro- 
static solution that fits the observed temperature estimates 
kT > 1.5 keV, if the only heat source is SN la events with 
a 'standard' rate. If the gas inside the cooling radius is in- 
deed so hot, there are two possible solutions of this problem: 
either repetitive gravitational heating is important due to 
violent changes in the gravitational potential of the cluster 
(or group) in a series of merging events, or the feedback star 
formation heated the gas up. To distinguish between these 
two or to find the proportions in which both are involved, 
the chemical evolution of the gaseous halo must be studied. 
Simple considerations for the gas enrichment in the course of 
hydrostatic recycling are presented in the next subsection. 

The recycling model represents a special exact steady- 
state solution of the general time-dependent hydrodynamic 
problem 



an equation for the ith element mass fraction Zi = Pi/p in 
the ISM: 



dp 

at 

d{pv) 



+ V • (pv) = ctp, - p ti , 

+ V(pv' 2 +p) = p\7(f> - ptiv, 



(24) 
(25) 



dt 

^+V.[v(E + p)] = 

ap*e, - p ti (E + p)/p - p 2 A + pv ■ V0 + V ■ q, (26) 

in which hydrostatic and thermal balance are decoupled. 
Here p = (7 — l)ep is the pressure, e is the specific internal 
energy, e* = c v Tb, the energy density E — p(e + v 2 /2), and 
q is the heat flux. Numerical solution of equations (|24|-(|2^) 
can be used to check the convective stability of the gas in 
the equilibrium recycling model and to study possible set- 
tling solutions in the inflow-outflow context. This requires 
multidimensional simulations and will be the subject of a 
separate paper. It is known from time-dependent numerical 
experiments on cooling flows that the distributed mass sink 
stabilizes the flow against the c ooling catastrohe for a suff i- 
ciently high condensation rate (Meiksin 1988; Friaga 1993). 
The general mechanism of stabilization and the domain of 
the parameter space where it operates have been described 
comprehensively in Section 3 above. Thus, numerical simula- 
tions of thermally unstable gas flows must take into account 
the condensational instability. This, however, mod ifies the 
gas dyn amics, because the gas becomes inelastic (Kritsuk 
1994 |). [ \ new type of elementary solution appears, repre- 
senting a sort of subsonic condensation wave, slowly propa- 
gating in the direction opposite to the local pressure gradient 
and 'reloading' the flow. The structure of wave fronts follows 
the random velocity field in the gas and therefore must be 
filamentary. Perhaps systems of bright filaments, observed 
in optical emission lines in some 'cooling flows', could be 
remnants of violent non-uniform condensation events in the 
core of the gas distribution during post-merging transitions 
to hydrostatic equilibrium. 



4.2 Coronal gas enrichment 

The hot gas metallicity is determined by the history of stellar 
mass loss and supernova activity in th e gala xy. Keeping the 
notation of Loewenstein & Mathews (1991), one can write 



dzi p, 

—r~ = — (a*u* 

dt P 



+ a sn y s 



(27) 



where y* ti and y sn ^ are the yields of the element, coming 
from stellar mass loss and SN la explosions, respectively. 
Under hydrostatic steady-state conditions equation (^) re- 
duces to 



a*y*,i + a sn y sn ,i 



y* 



~ysn,i 



and therefore the metallicity of the ISM may reflect that of 
the stellar population with the additional contribution from 
supernovae. The time-scale for gas enrichment, 



\ Zi J \ Zi 



(28) 



is proportional to the recycling time-scale t s , if all values 
in square brackets do not depend on the radius. Hence, if 
at some moment the hydrostatic gas with uniform Z has 
lower abundances than the stellar matter, a negative gradi- 
ent of metallicity is established in the course of evolution. 
It can be erased later, if the system has enough time to 
reach its 'chemical' equilibrium. This would not be the case 
if the stellar metallicity distribution (and, therefore, y*,i) in 
the galaxy were non-uniform. There are indications, based 
on magnesium indices as tracers of metallicity, that, while 
giant elliptical galaxies are inferred to have high metallici- 
ties, their metallicity gradients a re low, resulting in compa r- 
atively metal-rich stellar haloes (Thomsen & Baum 1989). 

The simplified picture above does not take into account 
some important effects, which can alter the metallicity dis- 
tribution and the absolute value of metallicity in the halo of 
a galaxy. First, if the enriched gas of stellar winds and su- 
pernova shells are not very well mixed into the hot IGM, 
but rather cool down after kinetic energy exchange with 
the hot medium, then there is no direct relation of the gas 
metallicity to the supernova rate or stellar metallicity of the 
galaxy. Secondly, the environmental effects could be impor- 
tant for gaseous haloes of dominant group/cluster members, 
as substructure merging in a galaxy distribution can heat 
the gas, and mix the enriched intragalactic medium with 
low- metallicity intragroup/cluster gas. While it seems to be 
possible to account for the first complication in the model 
proposed [just substituting a cooling function of higher Z in 
equation (hlh ] , one has to take care in applying the model 
to the hotter halo of a dominant member of a potentially 
dynamically young aggregation of galaxies, where the gas 
temperature still keeps a record of the recent merging event, 
but not one of the current SN activity. Adjusting the ther- 
mal equilibrium condition in this case would require an un- 
acceptably high supernova rate (and/or feedback heating) 
to get a higher value of To. 



5 DISCUSSION OF OBSERVATIONAL ISSUES 

It was pointed out in Section |^ that stellar surface brightness 
distributions, obtained with f3 = 1/2, are fitted well by the 
i? 1 ^ 4 -law. The scaling relation of the one-zone model showed 
that a mass deposition (star formation) rate proportional 
to p 2 allows stable recycling of the hot gas over the whole 
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galaxy with zero bulk motion of the gas. It will be illustrated ( [Thomas et al. 1986 ; Killeen & Bicknell 1988; Grillmair 



here th it the [3 = 1/2 model fits observed optical and X-ray 



surface brightness distributions well and reproduces major 
correlations. 



5.1 X-ray surface brightness distributions 

It is convenient to fit observed X-ray surface brightness dis- 
tributions E(i?) with a simple formula, which is a general- 
ization of the modified Hubble profile, 

21 -3/3 fit +0.5 



et al. 1994), the co oling radius is known to be ~ 68 kpc 
"(^"homas et al. 1986 ) , and the stellar ma ss estimate is M* (68 



kpc) ~ 10 12 M fljones fe Forman 1994| ). Therefore, accord- 
ing to equation (|23|), M(r) ~ 1.4 (r/68 kpc) ' 5 M yr _1 
for the adopted a = 4.7 x 1CT 20 s -1 . The mass deposition 
estimate in this case is close to values 2.2-3.0 M yr -1 , 



E(7?) = E 



1 + 



(!) 



(29) 



[see, e.g., Binney & Tremaine (1987)], where Eq, r c , and $ 



are free parameters and R is the distance from the centre 
in projection. Due to the different natures of X-ray ther- 
mal emission from diffuse gas and optical stellar light, the 
model with (3 = 1/2, which implies p» oc p 2 , provides optical 
and X-ray surface brightness distributions that match each 
otherj^. Therefore, /3fl t « /3fi t ,*- This is known to be the case, 

as far as relatively isolated elliptical galaxies are concerned heati ng: M s?s ~ 2.3 M yr _1 flMushotzkv fe Szymkowiak 



obtained by Thomas et al. (1986) for a cooling flow model. 
Similarly, for M87 the coo ling radius and mass estimates are 
nearly the same: 74 kpc (Stewart et al. 1984 ), and M*(70 
kpc) ~ 10 12 M ( Jones fe Forman 1994 ). The mean slope of 
the optical surface brightness distribution is —1.66, therefore 
M(r) ~ 1.4 (r/70 kpc) 34 M yr -1 . In this case the power 
index and the absolute value of the mass deposition rate 
disagree with the estimate, based on the 'standard' cool- 
ing flow analysis of Einstein observati ons: M(r) ~ 10 (r/7 
kpc) 2/3 M yr" 1 , for 5 <, r <, 70 kpc ( fttewart et al. 1984 - 
However, within the field of view of the Einstein Observa- 
tory SSS (~ 3 arcmin in radius) and FPCS (3 arcmin x30 
arcmin) our result M(3 arcmin ~ 13.1 kpc) ~ 0.8 M yr -1 
is in rough agreement with the upper estimates neglecting 



[see, foi example. Einstein IPC and HRI data on a sample 1986 ) and M, 



of six e larly-type galaxies, which includes three Virgo clus- 
ter mem bers: NGC 4472 (M49) NGC 4636, and NGC 4649 
(M60) ( |Trinchieri et al. 1986]); ROSAT PSPC observa- 



1.6 ± 0.3 M yr (Canizares et al 



1982), which have been accordingly rescaled. 



tions of NGC 4636 (flrinchieri et al. 1994|), NGC 43 65 and 
NGC 4382 (M85) flFabbiano, Kim fe Trinchieri 1994| ) ; anal- 



ysis of optical and X-ray surface photometr y of NGC 1399 



the centra l cD galaxy in the Fornax Cluster (Killeen & Bick- 
nell 1988|)]. 



In order to check if this is also true for M87, the X- 
ray dominant galaxy in the centre of a reasonably rich clus- 
ter, the results of 7?05^4T PSPC and HRI observations were 
analysed and compared with optical photometry of M87 at 
faint light levels in the B and V bands. The optical data 
are available to distances of ~ 260dis kpc F*| along the ma- 
jor axis and « 90 kpc on the minor axis ( [Carter fe Dixon 



5.2 X-ray gas temperatures 

When thermal instabilities occur in an accretion flow, two 
different spectral components are present in thermal X- 
ray emission, corresponding to terms pue (conventionally, 
softer X-rays from isobarically cooling comoving gas) and 
p 2 A (harder emission from the diffuse gas component) in 
equation ([io]). The spectral emissivity of such a composi- 
tion is 



e„(r) = p 2 (r)A„(T) + p ti (r)c v r„(T), 



1976 



The azimuthally averaged slope of the optical surface 
brightness profile is —1.66 and corresponds to (3m,* = 0.44, 
while /3flt for the averaged PSPC X-ray image is 0.45 ± 0.01 
(Bohringer 1994, private communication). Analysis of the 
NW sector of the HRI image, which is not contaminated by 
the radio components, gives (3m — 0.45 ±0.01. The compar- 
ison demonstrates a good agreement of X-ray and optical 
surface brightness distributions within the cooling radius, 
and similar orientations of the major axes of the outer opti- 
cal and X-ray isocontours of M87. 

The available surface brightness distributions allow one 
to determine the profiles of the mass deposition rate. Op- 
tical data give /3 flt = 0.42 for the cD galaxy NGC 1399 



II In this notation the surface brightness of stellar light follows 
equation ( p9| ) with /3fit,* instead of /3fj t , deprojected spatial dis- 
tributions of stellar and gas mass density follow p„ = p„ [1 + 

x 2]-3ft it- , _ x -60n%,* and p = Po [l + x 2 } -3/W2 ^ x -3/Sa t) re _ 

spectively, where x = r/r c , and r c could be different for the gas 
and the galaxy according to instrumental characteristics if the 
real core is unresolved. Note that f3 in ([]) and /3fl t in ( |2q| ) are 
different values. 

The distance of M87 is assumed to be 15 Mpc throughout this 
paper. 



where 



r„(T) 



7" 



1 



A„(T) 
A(T) 



dT 



(30) 



(31) 



( White & Sarazin 19871: ) . If the gas drops out due to ther- 
mal instability in equilibrium with stellar mass loss in the 
hot hydrostatic halo of a galaxy, the relative contributions 
from cooling diffuse gas and cooling condensates are fixed 
by the equilibrium conditions, and have similar spatial dis- 
tributions. In this case equation ( jjo| ) can be transformed 
into 

A 



e„(r) oc p*(r)A„(T) 



1 + 9 



TA, 



-P„(T) 



where for the adopted deposition law 
6(2 - u cq ) 



1 = 



6+1 



(32) 



(33) 



The radiative cooling function in equations (^) and (j3TJ) 
must take into account non-equilibrium ionization effects for 
temperatures less than several 10 6 K (see, e.g., Sutherland 
& Dopita 1993). 

While the idea of distributed mass deposi tion in cooling 
flows is widely exploited in physical models ( [Thomas 198(: ; 
White & Sarazin 1987a|; fearazin & Ashe 19891; iMeiksin 1986 ; 



Meiksin 199$ pavid et al. 1990| [Friaga 1993[ ), analysis of ob- 



served spectra is usually based on one- or two-temperature 
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fits, which appear to be satisfactory, mainly due to the lim- 
ited spectral resolution of the available detectors. This tech- 
nique does not account for the emission of cooling conden- 
sates, and temperature estimates, based on spectral infor- 
mation, suffer from limited spectral and spatial resolution 
of the detectors, model-dependent deprojection procedures, 
uncertain abundances (and their gradients) in the hot ISM, 
and unknown Hi column densities. 

Emission from condensates exceeds the diffuse emission 
at energies below ~ 2 keV when the background tempera- 
ture is 10 s K, and bel ow ~ 0.8 keV when the background 
temperature is 10 7 K (White & Sarazin 1987b). In the case 
of the ROSAT PSPC detector it will make a substantial 
contribution to the harder energy band of the detector, due 
to a blend of iron L-shell lines, if the gas condensations cool 
down from a temperature of about 3 — 4 x 10 7 K. In contrast, 
for background temperatures lower than ~ 10 7 K, the dom- 
inant contribution from isobaric cooling will be in the softer 



energy Dand. If spectral analysis does not take into account 



the emission from cooling condensations, the resulting pres- 
sure (and gravitating mass) estimates could contain errors 
of up to a factor of ~ 2 — 3. Consequently, the observational 
estimate of /3 spoc = cr* /(kT/fxm p ) is uncertain by the same 
factor. 

Reanalysis of the Einstein IPC, HRI, FPCS, and SSS 
observations of M87 has demonstrated the disagreement of 
single-phase models with mass determinations based on the 
optic al data, whi le the available X-ray data are explained 
well (Fsai 1994a). The total enclosed mass, derived from 
the model at r < 10 arcmin [i.e. right in the centr al 'cooling 
flow' r egio n , where the gas temperat ure drops (Bohringer 
et al. 1 )94; Nulsen & Bohringer 1995) to values lower than 
the kT = 2.05 ± 0.16 keV determined by Gmga (Koyama, 
Takano & Tawara 1991)], falls well below the mass estimate 
from the optical data. The multiphase model adequately 
explains the X-ray data and predi cts total ma sses, which 
agree with optical measurements (Tsai 1994b). The cen- 
tral gas temperature for this model is about 1.8 x 10 7 K. 
It seems that the main root of the problem is the lack of 
self-consistency in the standard procedure implemented for 
the gravitating mass determination in the 'cooling flow' re- 
gion, which is based on the hydrostatic assumption on the 
one hand and does not account for the distributed mass de- 
position (as the only means to keep the thermally unstable 
cooling medium hydrostatic) on the oth er hand. The best- 
fitting multiphase model of Tsai ( 1994b[ ) does not show any 
dramatic decrease of the temperature near the centre (which 
is typical for cooling flow models), but shows only a modest 
positive temperature gradient, dlogT/dlogr = 0.114, which 
is much smaller than that of the electron number density 
(dlogn e /dlogr = —0.491 in the core of radius r c = 6.64 
kpc and —1.36 at r > r c ). This can be considered as an ar- 
gument in favour of the isothermal approximation, adopted 
above, and justifies the usage of equation (^) for gravitating 
mass estimates. 



5.3 X-ray gas metallicity 

The distribution of iron abundance, as follows from the 
Ginga observation of the Virgo Cluster, gives a mean value 
for the central region within 1° of M87 of Z ~ 0.5 Zq. 
Because the angular profile of the abundance distribution 



is comparable to that of the Gmga collimator response 
(1°), the detailed structure remains unresolved. Broad Band 
X-Ray Telescope (BBXRT) observations of the interstellar 
medium in NGC 1399 indicated a mild positive tempera- 
ture gradient (kT = 1.0 — 1.1 keV for the central pixel of 4 
arcmin in diameter and kT — 1.1 — 1.2 keV for the outer pix- 
els at ~ 7 arcmin from the centre) and a subsolar metallicity 
with a negative gradient on the same scale, although the sta- 
t istical uncertainties do not preclude a constant metallicity 
( [Berlemitsos et al. 1993 ). 

Spatially resolved spectra, obtained with the ROSAT 
PSPC, can provide information on possible abundance gra- 
dients in the hot gas. Rapidly cooling regions in the centres 
of 'cooling flows' in NGC 1399 and NGC 4472 show an ap- 
parent trend of decreasing abundance with increasing radial 
distance (by a factor of ~2 for a scale of ~70 kpc), when 
T, Z, and the hydrogen column density, A?h, are tr eated as 
free parameters for the spectral fitting procedure (Forman 



et al. 1993; Jones fc Forman 1994). If JV H is held fixed, then 
the gradient is smaller ( |Forman et al. 199E ) . Similarly am- 



biguous results were obtained by Trinchieri et al. (1994) in 



their detailed study of ROSAT PSPC data on NGC 4636. 
If uniform 100 per cent cosmic abundances were assumed 
then two-temperature components were required to fit the 
observed spectra; if 20 per cent abundances were assumed 
then one component provided a fit of the same quality. New 
data from the Advanced Satellite for Cosmology and Astro- 
physics (ASCA) have confirmed the lower abundances for 
NGC 4636, however, with a mild negative gradient: from 
~0.3 at R a ~ 8.4 kpc to less than ~0.2 of the solar value at 
5R C (Mushotzky et al. 1994). Also there is good agreement 
between ROSAT and AS CA temperature profil es, T oc r~ - 2 
over a scale of ~ 30 kpc (Trinchieri et al. 1994). Even lower 
metal abundances, Z ~ 0.15 Zq, measured with ASCA were 
rep orted for NGC 1404 and NGC 4374 by Loewenstein et 
al. (|1994|). 



Although the observations do not give reliable informa- 
tion on the metallicity distribution in hot coronae, an as- 
sumption of uniform abundances with solar ratios, adopted 
above, could be an appropriate zero approximation for the 
region inside the cooling radius. Very low metal abun- 
dances detected by ASCA together with typical tempera- 
tures kT ~ 0.75 keV would imply incomplete mixing of gas 
injected by SNe. However, spectral fits based on the multi- 
phase gas model are needed to confirm the low abundance 
estimates. 



5.4 Lx—Lb correlation 

Using simple geometrical arguments, the scaling relation of 
the one-zone model (|l4|), and the virial condition po 



9ot/4:nG [see, e.g., Binney & Tremaine (1987)], one can 



write Lx oc poKr^ oc p, Cort oc p to o t r^ oc cr, r c . For galax- 



(Djorgovski & 



de Carvalho 199C), therefore Lx oc L^r c ' . Usi ng the 



Fabcr Jackson l aw for bright ellipticals [a, 



oc L° B 25 (|Fabe 



& Jackson 1976)] and the above relationship for the funda- 
mental plane, it is easy to eliminate r c and get Lx oc L^ 7 . 
Note that the Faber-Jackson relation can be derived from 
a simple self-regulated star formation scenario for a massive 



(M > 10 11 M ) spheroidal system (Lin fe Murray 1992) 



This scenario implies the initial expulsion of gas, combined 
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with that resulting from later stellar evolution, which may 
form a source of gas, feeding the hot haloes of individual 
galaxies. In agreement with this, observations show that the 
relation depends on the galaxy luminosity. While for the 
brightest (Mb ^ -20) ellipticals L B oc cr 4 - 2±0 ' 9 , for the less 
luminous ones (Mb <; —20) the deviations from the Faber 

<fT~ 

et al. Ifo82|). Thus ' — ' 1 "' < ' ' J 1 



Table X. Thermodynamic equilibria for hot coronal gas. 
Z= Z 



Jackso n la w may be considerable: Lb oc <j 2 ' 4±0 ' 9 (JUavies 
one has Lx oc L s 65 and Lx oc L s ' lb for 
bright and faint ellipticals, respectively. These results can 
be compared with Lx oc L^ 7±0 ' 3 (Canizares, Fabbiano & 
Trinchieri 1987) and Lx oc l| 18±0 ' 20 (Donnelly, Faber & 
O'Connell 1990). Note that the derived relation of X-ray to 
optical luminosities is quite different from the naive scaling 
Lx oc Lb for SN- heating dominated models. 

In the framework of the recycling model the deviations 
from the 'mean' X-ray to optical luminosity relatio nship [the 
observed rms residuals of logLx are 0.33 — 0.35 (Donn 
et al. 19^9C|)] may be attributed to the combined effects of: (i) 
galaxy-to-galaxy differences in the spectra of perturbations, 
which stimulate thermal instabilities in the hot gas; (ii) dif- 
ferent rates of SN activity and/or stellar mass loss rates; (iii) 
variations of the chemical composition of the hot gas. Each 
of these circumstances determines the equilibrium gas con- 
tent of a galaxy and most of them may be environmentally 
influenced (see also Section 5.6 below). The most dramatic 
deviations may occur if the efficiency of condensation b be- 
comes lower than its critical value 6 cr it • This could s timulate 
oscil lati ons of the gas density with amplitude ~5 (Kritsuk 
1992|)~^nd, therefore, initiate the random walk of a galaxy 
about the mean value of Lx- The possibility of such desta- 
bilization is one of the natural explanations for considerable 
scatter in the logLx — log Lb plane. 



5.5 Lx— T correlation 

For the /3-model T oc a 2 and therefore, using the Lx — Lb 
relation above and the Faber-Jackson law, one obtains Lx oc 
T 3,4 . For the high- and low-luminosity ellipticals Lx oc y 3 - 44 
and Lx oc T 2 54 , respectively. Note that this scaling relation 
is similar to Lx oc T 3 3 , for a sa mple of galaxy clusters ob- 
served by Ginga (Arnaud 1994). For cluster cooling flows 
the observed X-ray luminosity-temperature relation has the 
same scaling, but depends also on the mass deposition rate 
(estimated with the 'st andard' cooling flow analysis tech- 
nique): L x oc T 3 ' 3 M°- 4 (|Fabian et al. 1994|). 



5.6 Gas content and environmental effects 

The relation of the SN-associated heat source temperature 
To to the equilibrium gas temperature and density, as well 
as Z and b,is illustrated here in Table 1, computed using 
equations (^), ( |l6| ) a nd (|j ^) and the cooling functions of 
Sutherland & Dopita (1993). It shows that higher metallici- 



ties result in lower To values and lower gas densities (X-ray 
luminosities) for a fixed temperature of the hot diffuse com- 
ponent of the ISM. Also the higher condensation efficiency 
provides a lower Lx- When the efficient feedback star for- 
mation is triggered by the recent merging event, a higher 
equilibrium X-ray luminosity is expected. This might be the 
generic case for dominant galaxies in rich and, therefore, 
more active environments. In relatively isolated galaxies, 
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The heat source temperature Tq 7 and the equilibrium gas tem- 
perature T 7 are given in units of 10 7 K. The integration constant 
C—15 in equation (rj) is computed for the stellar mass loss rate 
a = 4.7 X 10 -20 s _1 and given in units of 10 -15 (g cm" 3 )!. 



as is illustrated in the following, the standard value of the 
SN rate can readily explain the observed gas temperatures. 
Note, however, that the uncertainties of the To estimates 
inherit the uncertainties of the model cooling functions and 
their first derivatives. Therefore the values of To may be only 
30-50 per cent accurate. 

Recent ASCA observations of bright elliptical galax- 
ies in the Virgo Cluster revealed hard X-ray components 
with colour temperature kT > 2 keV in addition to ex- 
t ended thermal X-ray em ission of temperature kT ~ 1 keV 
(Matsushita et al. 1994). This hard component is primar- 
ily attributed to the integrated emission from discrete X-ray 
sources, associated with the optical galaxy. For galaxies lo- 
cated in the area of extended cluster X-ray emission, the dis- 
crete sources alone cannot account for all the hard X-ray flux 
detected with ASCA. It was assumed, therefore, that the 
hard components of NGC 4406 (M86) and NGC 4374 (M84) 
are contaminated by the foreground/background Virgo in- 
tracluster medium emission. In the case of NGC 4472 the 
hard-band image is elongated and displaced from the op- 
tical galaxy and from the soft-band image as well. This is 
interpreted as the locally enhanced ICM emission plus the 
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discrete-source contribution, although their respective con- 
tributions remain unspecified. The temperature of the soft 
emission from NGC 4472, kT = 0.83±0.04 keV, corresponds 
to an SN heating temperature To w 3.4 x 10 7 K for the ob- 
served metallicity Z ~ 0.42 Zq and the condensation effi- 
ciency b = 0.7. For NGC 4636 located at the outskirts of 
the Virgo Cluster, ROSAT PSPC data imply a similar hard 
(kT > 0.9 keV) component, which is more extended than 
the soft one (kT ~ 0.6 — 0.8 keV) and probably is not di- 
rectly related to the galaxy. However, this is only true if 
100 per cent cosmic abundances are assumed. With 20 per 
cent cosmic abundances only one temperature c omponent 
(kT ~ 0.6 — 1.0 keV) is required for spectral fits (Trinchieri 



et al. l|9J). 

The BBXRT d&t& for NGC 1399, a cD galaxy in a poor 
cluster, suggest a gas temperature kT = 1.0 — 1.2 keV and 
metalli cities of abou t 0.5 Zq with large uncertainties (Ser- 
lemitso > et al. 1993). The analysis of two-temperature mod- 



els revealed (similarly to the Virgo galaxies) a hard compo- 
nent of X-ray emission with kT ~ 6 keV and provided a 
lower temperature ~ 1 keV and a higher metallicity 0.76 Zq 
for the primary soft component. This would imply a char- 
acteristic heat source temperature To ~ 4 x 10 7 K. 

An accurate deprojection analysis of the ROSAT PSPC 
spectral data on M87 provided temperature estimates for 
the centre of the 'cooling fl ow' region (r < 3.33 arcmi n ~ 15 
kpc) offcT = 1.1-1.4 keV flNulsen fc Bohringer 1995| ). Such 
a temperature would require To ~ 5 x 10 K for the subso- 
lar metallicity Z ~ 0.5 Zq and an even lower value for solar 
metallicity. The condensation efficiency b = 0.7 roughly cor- 
responds to q = 0.874 - the value for the best-fitting mul- 
tiphase model of Tsai ( 1994b ), see equation (^j|). Note that 
t he te mperature estimate was made in Nulsen & Bohringer 
(1995) for the region where the errors due to contamina- 
tion by the ICM emission are small. The straightforward 
approach would require To ~ 9 x 10 7 K for the gas temper- 
atures detec ted with wide-beam instruments [e.g. kT ~ 2 
keV, Ginga (Koyama et al. 1991)]. Note, however, that the 
multiphase model based on Einstein data gives a central 
temperature estimate of T ~ 1.8 x 10 7 K. 

Although in all three cases discussed above the low gas 
temperature estimates were obtained on the basis of a single- 
phase description, if they are considered as rough estimates 
for the temperature of the diffuse component of the cool- 
ing gas, the standard SN la rate appears to be sufficient 
for thermodynamic equilibrium at the measured temper- 
atures. The analysis of ASCA observations, based on the 
multiphase model, would help to check whether the diffuse 
component of the multiphase ISM in central dominant galax- 
ies actually has a low temperature ^ 1 keV and nearly solar 
metallicity. In this case a kT ^ 1 keV galactic atmosphere 
might be embedded in the extended intracluster gas with 
a higher temperature and lower metallicity. ASC A data on 



NGC 4696/Centaurus lend support to this option (Fukazawa 



et al. 1994). Such a complex quasi- hydrostatic configuration 
can be described with a hybrid model, including the recy- 
cling interior of the dominant galaxy matched at the cooling 
radius [where t c (r C ooi) — Hq ] with an isothermal /3-model 
for the hotter tenuous intracluster gas. If the ratio of the 
stellar-to-gas mass is low and the cooling radius is larger 
than the extent of the optical galaxy, then a cooling flow 
is unavoidable. It starts from the cooling radius, enters the 



galaxy, and disappears, providing material for star forma- 
tion. 



6 CONCLUSIONS 

A purely hydrostatic model for the hot gas in normal el- 
liptical galaxies has been presented as an alternative to the 
'standard' cooling flow picture. The model ind epend ently 
develops an option suggested by Thomas et al. (1986) that 
perhaps there are no radial gas flows, so that all gas present 
at some radius originated there (e.g. from stellar mass loss) 
and cools and is deposited in situ. Energy injection from SNe 
and feedback heating due to star formation from deposited 
material have been also incorporated in the description of 
equilibrium gas recycling. Formally, the model is quite simi- 
lar to the st eady-state cooling flow formulation of Sarazin & 
Ashe ( 1989 ) but with zero bulk motion of the gas. (Note that 
the notion of the cooling flow itself implies highly subsonic 
velocities and there are no observational data that would 
directly confirm the existence of the flow.) Therefore, the 
recycling model, representing a 'minimal cooling flow', gives 
the exact lower estimates of the mass deposition rate in el- 
liptical galaxies and fills the gap in a family of models for 
hot gaseous coronae, being a fully self-consistent hydrostatic 
model which for the first time takes into account the whole 
set of sources and sinks of mass and energy. It has the ad- 
vantage of an analytical description of gas equilibria, which 
allows the study of their linear stability. However, the global 
thermal-convective stability of such a hydrostatic recycling 
atmosphere of a galaxy must be the subject of a separate 
numerical investigation. 

Summarizing the results of this study, it is worthwhile 
to stress that the galaxy, located in the centre of the cooling 
region, may play an important role, providing a complex 
mechanism for hydrostatic support of thermally unstable 
gas. The moderate amount of SN la heating in combination 
with the distributed mass sink is able to compensate the ra- 
diative cooling and the stellar mass loss in such a way that 
quasi-isothermal hot gas is kept in both hydrostatic and ther- 
modynamic equilibrium. The latter becomes possible due to 
a first-order phase transition, which separates the stable hot 
gas phase for a sufficiently high efficiency of condensation. 
Having no energy input and/or mass sink it is impossible 
to avoid a cooling catastrophe in the centre of the 'cooling 
flow' region, where radiative cooling and stellar mass loss 
are efficient. 
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